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Supplementary Material 
The melting profile 

To measure the melting curves of DNA by spatially 
scanning the thermal profile, we assumed that the DNA 
was at thermal equilibrium at each point of space. This 
assumption is true if thermal equilibration is fast com- 
pared to the drift and diffusion velocity. The drift can 
originate from the electrophoretic mobility of the DNA 
or its convection in a electroosmotic flow (if present). 
Because we use AC voltages the average value of the 
drift is zero. Notice that most of the melting profile 
is measured outside the hole region where the elec- 
tric field is small; only the small ^-values correspond 
to the vicinity of the hole where the field is apprecia- 
ble. The excursion of the molecule in the temperature 
profile during one period can be evaluated as follows. 
The electrophoretic mobility of DNA is largely inde- 
pendent on the DNA size. Its value 1 is 3 x 10~ 4 cm 
V -1 s _1 . The maximum electric field can be estimated 
by dividing the maximum voltage by the characteristic 
size, that is E = 25/50 x 1(T 6 = 5 x 10 5 V/m. The elec- 
trophoretic velocity of the DNA in the hole is then equal 
to 1.5 x 10~ 2 ms -1 . Hence the drift distance during half 
a period of the field is of order 1 fim. At larger distances 
from the pore the electric field decreases as 1/r 2 and the 
drift distance is correspondingly decreased. 

We could not detect any electroosmotic flow. If any 
such drift existed it would have been localized at the 
hole where all flow lines converge. We used a Teflon 
septum in a 1M KC1 solution. This means that the 
zeta potential and therefore the electroosmotic flow 
are largely reduced. The electroosmotic coefficient^ for 
PTFE (Teflon) capillary tubing is below 1(T 4 V^ 1 s" 1 
rendering this drift negligible. 

Lastly we point out that the temperature profile is 
also estimated using DNA molecules with similar elec- 
trophoretic mobility and undergoing the same electroos- 
motic flow. Hence, the temperature that we measure at 
a fixed point is also the mean temperature over the ex- 
cursion of the DNA in the temperature profile. 

Applicability of the device 

This device can be a useful tool to produce large tem- 
perature gradients. Several parameters may affect the 
applicability of the concept: 

1 - AC frequency. Heating occurs at all frequencies 
and increases with decreasing frequency. At small fre- 
quencies, however, the drift of the DNA molecules in 
one period is diminished. We tested frequencies be- 
tween 1 to 30 kHz without any significant change in 
the results. At large enough frequencies the current 
across the membrane can be dominated by the capaci- 
tive charging of the membrane. The crossover frequency 
is oj = (RhCs)- 1 = arl/{aC s )- x , where R h ~ 2300 Q 
and ro are the pore's resistance and radius, C s ~ 0.4 pF 
is the septum capacitance, a is the electrical conductiv- 
ity and a is a geometry-dependent factor. For the 50 
/jm thick membrane we used the crossover frequency is 
1 GHz. For thin membranes capacitive charging should 



be considered and in this case the ion flux occurs every- 
where in the solution and the heating is consequently 
less localized and more homogeneous. Note however 
that the product RhC s does not dependent on the mem- 
brane thickness for a planar capacitor and a cylindrical 
hole. 

2- Buffer composition. In the Teflon membranes we 
used the values of pH and salt concentration do not dra- 
matically affect the zeta potential at the pore surface. 
Other polymers or solid state membrane can be affected 
by this parameters. 

As stated above, heating occurs within a few tens of 
microns of the pore and the Peclet number is small. 
In this limit, the analysis and heating profile should 
be unchanged even in the presence of convection. No 
changes were detected in our experiments with solutions 
of concentrations in the range from 0.2 M to 1M. 

Numerical calculations 

The governing equation for the salt concentration 
C^, electric potential (j>, and temperature T 



- V(L' ± VC ± ± ez^C^cj)) 
<9 t (V(eV0) + V(ctV0) 
pC p d t T - V(kVT) - er(V(/>) 2 



(1) 
(2) 
(3) 



were solved using a finite element method with a high 
resolution mesh near the pore. C ± , Z? ± and /i ± are 
the concentration, diffusion coefficient and mobility of 
positive ("+", K) and negative ions ("-", CI) in water, 
respectively. The diffusion of K + and Cl~ ions were 
taken to be D+ = 1.9 x 10~ 9 m 2 /s and D^2x 10~ 9 
m 2 /s and the mobilities depended on the diffusion via 
Einstein's relation, e is the electron's charge and the 
valency z was chosen to be z — 1. The first equation 
was calculated in the aqueous region only as ions can- 
not penetrate the hard walls. The second equation was 
solved in the aqueous regions with electrical conduc- 
tivity given by o — e 2 (p + C + + p~C~) and dielectric 
constant s = 80eo, where Eq is the vaccum permittivity, 
and with a = and e — 2.\eq in the solid parts. The 
usual boundary conditions at the interfaces between di- 
electric or conductive materials were used. 

Due to the total electro-neutrality and the AC fields 
we used, convection can be neglected and the Navier- 
Stokes equation need not be solved. Subsequently the 
heat equation (third equation) has only Joule heating 
terms and no viscous dissipation. C p , k and p are the 
specific heat, heat conductivity, and mass density of 
the aqueous region and Teflon cell. We used C p = 4180 
J/Kg K°, k = 0.58 W/mK°, and p = 1000 Kg/m 3 
in water and C p = 1000 J/Kg K°, k = 0.27 W/mK°, 
and p = 2140 Kg/m 3 in Teflon. At the water- Teflon 
interface we conserved the heat flux, and at the edges of 
the cell we used Newtonian heat loss to the environment 
(proportional to T — T^,) with temperature at infinity 
Tqo equal to the room temperature. 

iE.Stellwagen, N.C.Stellwagen, Electrophoresis 23, 2794 (2002). 
2 P.Camillery, Capillary electrophoresis, Ed: CRC Press LLC), 
103 (1998). 
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Localized Joule heating produced by ion current focusing through 
micron-size holes 
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We provide an experimental demonstration that the focusing of ionic currents in a micron size hole con- 
necting two chambers can produce local temperature increases of up to 100° C with gradients as large as 
1° K/j,to _1 . We find a good agreement between the measured temperature profiles and a finite elements- 
based numerical calculation. We show how the thermal gradients can be used to measure the full melting 
profile of DNA duplexes within a region of 40 /im. The possibility to produce even larger gradients using 
sub-micron pores is discussed. 
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The creation of local heat sources and large thermal 
gradients in confined aqueous environments is a chal- 
lenging problem due to the rapid diffusion of heat in 
water. Several solutions were proposed, such as heat- 
ing micro/nanoparticles^ by magnetic induction or by 
using a focused laser beam 3 -^. Local thermal gradi- 
ents in microchannels were used to sort and concentrate 
molecules. Several approaches were developed such 
as thermophoresis^, Temperature Gradient Focusing^, 
Field Gradient Focusing^, and isoelectric Focusing^. 
These techniques use high DC voltages (tens to hun- 
dreds of Volts) and thermal gradients in the range of 
0.01° K/xm - . The typical volumes are micro to nano- 
liters. In this letter, we show that a few tens of picol- 
iters can be strongly heated by focusing an ionic current 
through a micron size hole in a saline solution. The re- 
sulting gradients are of the order of 1° K/zm _1 . 

We use a custom made cell composed of two cham- 
bers separated by a 50 /im thick Teflon septum (see 
Figure [T]). A conical hole (40° half angle, minimum 
diameter ro = 7.5 /im) is punctured in the center of 
the septum. The AC voltage (~ 100V PP , 10 kHz) is 
applied across the chambers using platinum electrodes. 
We use tris buffer, 1M KC1, pH 7.4, of electrical conduc- 
tance a = 107 mS/cm. The DC value of the electrical 
resistance Rh across the hole is 2300 f2. The septum 
capacitance is estimated to C s = 0.4 pF. It results that 
the current across the chambers is mostly resistive at 
10 kHz. The Joule heating power Pj dissipated in the 
hole is proportional to the mean root square current 
i rms . We measure the local temperature profiles along 
the vertical hole axis in the lower chamber. The tem- 
perature is derived from the confocal detection (A = 532 
nm) of the calibrated fluorescence of TetraMethylRho- 
damine (TMR) grafted at the 5'end of DNA oligomers. 




FIG. 1. A: Schematic representation of the experimental 
setup. The two chambers are separated by a 50 fim thick 
Teflon septum punctured by a 7.5 /im radius conical hole 
B: Close up on the hole region. C: Model set up used for 
our finite element modeling. D: close up on the pore region. 
The ion current lines are represented and the map of the 
temperature profile is color coded. 



Large enough currents, i rms > 1.9 mA, result buffer 
vaporization. Smaller values of « rms lead to station- 
ary thermal profiles within a few seconds. Figure [2] A 
shows the temperature increase AT = T — T^, where 
Too = 298°K is the room temperature, along the verti- 
cal axis for i rms = 1.81 mA. Over the first 40 /im , the 
average temperature gradient is 1° K/im -1 . As a first 
approximation we model the heating power P e g as be- 
ing distributed uniformly in an effective sphere of radius 
r c ff. The temperature then reads: 



AT(r) 
AT(r) = 
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a ' Electronic mail: virgile.viasnoff@espci.fr 
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where r is the distance from the hole center and k = 0.6 
Wm _1 K _1 is the thermal conductivity of water. Fig- 
ure OA shows the measured temperature profile for hole 
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radius tq = 7.5 /im and power Pj = 6.5 x 10 -3 W. The 
best fit parameters for our model are r c g ~ f7.8 /im 
and a P e ff = 5.2 x 10~ 3 W. This model predicts that 
when r > r c s the value of (r/ro)AT(z)/AT(ro) is a 
constant independent of ro and Pj . Figure shows 
that the curves for several values of ro and Pj can be 
scaled provided that we take r c ff = 2.4rg. 
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FIG. 2. A: Profile of the temperature increase AT(z) along 
the vertical axis below the hole center. The pore's radius is 
7.5 (jm and the current is i rms = 1.81 mA. The average ther- 
mal gradient over the first 40 fim is 1° K/im^ 1 . At distances 
larger than the hole radius the temperature profiles AT(z) 
scale as 1/z. The dotted line is the fit obtained when mod- 
eling the heating source by a spherical source of uniform 
heating power Pj = 5.2 x 10 -3 W and radius r c g — 17.8 
/im. The solid line is the finite element calculation. B: 
Rescaling of all temperature profiles obtained with a hole of 
7.5^tm (open circles) and 20/im (full triangles) for various 
currents (0.5mA < i rms < 1.85mA). We used r e ff = 2.4ro. 
The straight line corresponds to the scaling of the analytical 
model. Inset: Values of AT(r e ff) as a function of I rm s for 
the hole of 7.5/mi (open circles) and 20pm (full triangles). 



In order to explore the influence of the pore geometry 
on thermal profiles we use a finite elements approach 
with various hole radiuses ro and pore lengths L . The 
numerical calculation accounts for mass transport ([2]), 
electrostatic potential (|3|) and heat dissipation (|4]) in 
the following way: 

5 f C ± - V(L» ± VC ± ± ez/i ± C ± V^>) = , (2) 



d t {V{eV(f>) + \7{<j\7<f>) =0 , (3) 
P C p dtT - V(kVT) - d(V0) 2 = , (4) 

C ± , D ± and fj,^ are the number densities, the diffu- 
sion coefficients and the electrophoretic mobilities of 
the positive and negative ions, respectively. <fi is t ne 
electric potential, e is the local dielectric constant, a 
is the electrical conductivity, e is the electron's charge, 
and C p is the heat capacity. The steady-state solutions 
were obtained with an AC voltage at the electrodes. See 
Supplementary Material for details. 
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FIG. 3. Steady state temperature gradients and profiles 
(inset) below the narrowest constriction of the hole given by 
finite elements calculation. V z T(z) and AT(z) are scaled by 
the average heating power Pj for various ratios of hole size 
ro = 10 /im with the membrane thickness L. The values of 
r /L are 0.01 (o), 0.05 (a), 0.1 (□), 0.2 (o), 0.5 (+). AT{z) 
agrees with Eq. [1] and decays like z~ x at z > ro for ro/L 
larger than 0.1. For lower ratios the thermal gradients are 
less steep. 



Figure [3] shows the calculated temperature AT(z) 
and gradient profiles V z T(z) along the z axis at a fixed 
heating power Pj for several values of the aspect ratio 
ro/L. When tq/L > 0.1 the temperature is well de- 
scribed by AT(z) r~j z^ 1 for z values larger than 2ro- 
In the limit of small aspect ratio the heating is mainly 
localized in the hole and smoothly spreads in the lower 
chamber producing smaller temperature gradients. We 
conclude that ro/L should ideally be between 0.1 and 
1 to obtain the sharpest gradients and largest temper- 
ature increases. 

The experimental thermal gradients generated by a 
hole of diameter ro = 20 /jm are used to determine 
the melting profile a DNA duplex: strand 1: TMR 
5'-TCAGACCG(TC)i 5 -3', strand 2: 5'-CGGTCTGA- 
3' IowaBlack. The DNA was gel purified to obtain a 
95% hybridization efficiency. The fluorescence inten- 
sity of the TMR is quenched 20-fold on average upon 
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hybridization with IowaBlack. With a proper baseline 
calibration the fluorescence intensity measured at the 
laser spot can be used to determine the fraction of hy- 
bridized duplexes. We measure the fluorescence pro- 
files at various values of z rms for strand 1 only, and for 
the hybridized duplex with a 1:1 ratio of both strands 
(see Figured]). Assuming (i) a local thermal equilib- 
rium, (ii) a two state model where the 8 mers are ei- 
ther fully hybridized or completely open, and (iii) an 
efficient quenching for all temperatures, we can extract 
the dissociation coefficient a of the hairpin as a func- 
tion of temperature and distance from the pore. Fig- 
ure @] shows the full melting profile obtained over a dis- 
tance of 45 /im for a hole of radius ro = 20 /itm. Fol- 
lowing the typical melting curve analysis for bimolec- 
ular equilibrium^, we extract the thermodynamical 
melting parameters of the DNA structure. We find 
AH = -66 ± 6 kCal mole" 1 and AS = 182 ± 25 
Cal mole -1 K -1 , in good agreement with the thermo- 
dynamical parameters calculated with MFold 1 "^ under 
similar salt conditions (AH = —62 kCal mole -1 and 
AS = 169 Cal mole -1 K -1 ). Spatial temperature gradi- 
ents have the advantage over traditional melting curves 
techniques that all temperatures can be probed simul- 
taneously. This method works if the diffusion and/or 
drift of the DNA molecule across the thermal profile 
is slow enough to allow thermal equilibration. As de- 
rived in the Supplementary Material, electrophoretic 
and electroosmotic drifts are negligible in our exper- 
iments. The Brownian diffusion coefficient 1 ^ for our 
DNA molecules is of order Z?dna — 10 -7 cm 2 s -1 . The 
distance over which the temperature changes by 1°C 
is d=l /im. The diffusion time across this distance is 
t = gP/Ddna = 0.01 s. Since small molecular beacons 
are reported to open over a characteristic time of 10 -4 
s 1 ^, the approximation of local equilibrium is satisfied. 



the limit of small Peclet numbers i.e if the convective 
velocity is small: v < n/pC p ro = 10 -5 /ro m s -1 . In 
this limit thermal quenching rate of 1° K//is can still 




40 50 60 70 80 

T / °C 



FIG. 4. The dissociation ratio a = 

-Zfl U o(DNA + quencher)/7fl uo (DNA only) as a function 
of distance and temperature for i rms — 2.9 mA. A decrease 
of 90% is obtained over 45 fim or 35° C. The hole diameter 
is 40 /J.m. Inset: fluorescence intensity as a function of the 
ionic current for a single DNA strand labeled labeled with 
TMR only (full symbol) and a DNA duplex with TMR 
quenched with IowaBlack (open symbol). The lines are 
guides to the eyes. 



be achieved. See Supplementary Material for a detailed 
discussion. We believe that this approach facilitates 
the creation of large thermal gradients in sub-micron 
regions with potential applications for fast denatura- 
tion and thermal quenching 1 ^ to study local chemical 
reactions. 



We briefly discuss the influence of various parame- 
ters on the applicability of our device. For frequencies 
w < (RhCs)^ 1 = ar^KaCs), where a is a geometry- 
dependent factor, the ionic current is mostly resistive. 
It is reasonable to assume that electroosmotic flow are 
negligible due to the use of AC fields. Electroosmosis 
can be significantly enhanced in smaller (~ 100 nm) 
or charged pores when current rectification occurs 1 ^. 
In the lower chamber natural convection is minimized 
since the hot spot (the hole) is located above the cold 
region. In addition, there is negligible fluid transfer 
between the upper and lower chambers. The small ex- 
tent of the lower chamber also increases the instability 
threshold for natural convection. This situation con- 
trasts with focused laser heating for which convective 
rolls are observed along the optical axisi Convection 
can be externally applied to drive the molecule through 
the thermal profile. Our conclusions would still hold in 
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